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Abstract 

The colour connection structure of QCD (2 — > 2) processes is discussed, 
with emphasis on its apphcation to the supersymmetric 2 parton 
2 sparton processes, which are currently being implemented in the 
HERWIG Monte Carlo event generator. The procedure described by 
Marchesini and Webber is found to be inadequate, and a new method 
is proposed. However, this alteration is unlikely to significantly af- 
fect the theoretical predictions for soft gluon radiation. A complete 
list of supersymmetric QcD 2 — > 2 matrix elements and their colour 
decompositions is presented. 
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1 Introduction 



The simulation of soft gluon radiation in hard (supersymmetric) QCD pro- 
cesses |jl|, H, ^ H requires that the corresponding matrix elements be rear- 
ranged according to the colour connections (defined by the colour flows) in 
the process [IH . In brief, this is because the colour connections in the parent 
process determine the cones in which the soft gluons radiate and hence the 
hadronisation occurs. 

This rearrangement of terms is automatic if there is a unique colour flow 
associated with the process, as is the case for the QcD qq' — > qq' scattering 
whose colour flow is shown in figure |I], but for more complex processes the 
procedure involves some ambiguity. The purpose of this paper is to discuss 
and analyse this ambiguity, to propose a consistent and practicable method, 
and to illustrate the application of this new method in supersymmetric QcD 
2 parton 2 sparton processes which are being implemented in the HERWIG 
Monte Carlo event generator 0, |^, |^ . 

1.1 Colour flows in QCD 

The matrix elements for processes with more than one colour flow consist of 
the 'planar' terms and the 'nonplanar' terms. The planar terms are those 
with single colour flows and the nonplanar terms are those with no single 
colour flow. The nonplanar terms are always suppressed by some inverse 
powers of Nc- 

The colour flows for the four distinct 2^2 QcD processes are shown in 
figures For concreteness, let us consider the QcD process qq — > gg, for 
which the leading-order spin- and colour-averaged matrix element squared is 
given by 0: 

The first term in braces is a planar term corresponding to the t-channel 
colour flow, and the second term corresponds to the w-channel colour flow, 
as depicted in figure ^ The third term, suppressed by the factor (l/N^), is 
a nonplanar contribution corresponding to a mixed colour flow, which needs 
more care in its treatment when we consider the rearrangement according to 
the colour connection. 

Note that apart from the overall gauge invariance and positive definiteness 
of the matrix element squared, each of the above three terms is also gauge 
invariant and, in the case of the planar terms, positive definite. The gauge 
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invariance follows from the fact that colour is formally an observable at Born 
approximation. As for the positive definiteness of the planar terms, this is 
obvious since the modulus squared of any part of the full matrix element 
must also be positive. This gauge invariance allows us to uniquely identify 
the planar terms and the nonplanar part in each process. 

1.2 Radiation from colour connected partons 

Let us recall the results of Ellis, Marchesini and Webber [||, § concerning the 
coherence of soft gluon radiations [^, ^ in hard QCD processes. 

For each colour flow (each planar term) there is a cone around each in- 
coming and outgoing parton direction bounded by the angle between the 
parton and the parton which is colour connected to it[|. In the case of gluons 
(and gluinos) there are two such cones, one for the colour and one for the 
anticolour. The cones define the bounds for the soft gluons to be radiated (in 
the angular ordering approximation) and hence the hadronisationQ to occur, 
to leading order in Nc- 

As for the nonplanar part, this can be distributed in any way between 
the colour flows. Let us consider this in detail. 

Returning to the example of qiq2 — > q^Qa introduced earlier, the radiation 
pattern, in the notation of is: 

(2) 
with: 

Wt = Ca[Wm + Wi^ + W2a-Wi2]+2CfWx2 (3) 
Wu = Ca[Wsa + Wu + W23-Wi2]+2CfWu (4) 

W^n.p. = CA[Wi3 + W2i + Wu + W23-2Wi2]+2CFWi2. (5) 

Cp and Ca usual the QcD colour factors (A^^ — l)/2Nc and (Nc) 

respectively. Terms proportional to Ca correspond to the radiation from 
the gluon legs, whereas the terms proportional to 2Cf correspond to the 

^This statement has not been explicitly verified in the supersymmetric case involv- 
ing either long-lived spartons in the final state (the light gluino, whose existence is still 
controversial — see |l^.) or short-lived ones in the intermediate state. 

^When the (s)partons are massive, the analysis needs some modification [0. For 
the heavy quark case, there will be a 'screening' of the coUinear direction For the 
supersymmetric case, the analysis has not yet been carried out. 
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radiation from the quark legs. Wij are the radiation functions for the parton 
pairs {i,j} and are defined as: 

Wi ■ = " (6) 

for the emission of a soft gluon of momentum pg. In terms of the angles 
between parton pairs, 



W-=—-^f- (7) 

^g ^"-9 



where = 1 — cos 9ij = Pi ■ pj/EiEj. Eg is the soft gluon energy. This 
expression contains poles at 6ig = and 9jg. These correspond to the two 
collinear directions in gluon emission. Explicitly, 

W., = W-, + W-^ = ^{^ + f | + (8) 

W^j have collinear singularities only at 6ig = 0. It can be shown that the 
azimuthal average of W-j around the i-th parton direction is a step function 
with cut-off at 6ij, hence the above claim concerning the radiation pattern. 

We now consider modifying the planar terms such that the sum of them 
is equal to the original matrix element squared. Note that: 

Wt-W^,p, = Ca[Wu + Wu-Wu-W23] (9) 
W^„-W^n.p. = Ca[Wu + Wu-W,s-W24]. (10) 

In both expressions the radiation cancels in all four collinear directions, and 
so it is reasonable to approximate W^.p. by either Wt or Wu- 

Hence it can be deduced that the radiation due to the (1/iV^) suppressed 
nonplanar term can be treated approximately by distributing this term be- 
tween the two colour flows in some ratio. 

From the viewpoint of practicality in Monte Carlo simulations, this dis- 
tribution should be such that the modified planar terms, which are called 
'full terms' in should be positive definite. 



2 The MW procedure 

The procedure of Marchesini and Webber (MW) for evaluating the full terms 
is as follows |]l|: 

• the full term should have the same pole structure and crossing symme- 
try as the planar term; 
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• it should remain positive definite in order to be interpreted as a prob- 
ability distribution; 

• the sum of the full terms should give the exact lowest order 2 — > 2 
matrix element squared. 

The second and the third criteria above are essential as argued earlier. 
As for the pole structure, the introduction of an extra physical pole is bound 
to drive at least one of the full terms negative for some permutation of the 
external legs so this is a natural consequence of the second criterion. The 
requirement of correct crossing symmetry is never utilised in practice, and it 
is in fact too constraining to be practicable, as will be shown later. Discarding 
this requirement, we have essentially only the positivity and the sum to 
constrain the colour rearrangement, which is not sufficient to determine it 
uniquely. One might therefore introduce the following criterion, which is 
similar to the requirement of (not only physical) pole structure, and can be 
regarded as being inherent in the MW procedure: 

• when a nonplanar term contains poles corresponding to two colour 
flows, this is split up by partial fractions, viz 

111 

- + - + - = 11 

st su tu 

for massless partons. 
There are several disadvantages associated with this method: 

• the decomposition of the nonplanar part is not unique without the 
additional criterion introduced above, and with this criterion the pro- 
cedure sometimes fails and/or is still not unique, depending on the 
exact procedure by which the partial fractions are split up; 

• the decomposition is not general under permutations of external legs; 

• the decomposition becomes laborious when the number of colour flows 
is increased, and verification becomes an impossible task. The number 
of colour fiows is a factorial function of the number of external partons. 

The first two of these points are illustrated well for the case of the supersym- 
metric QCD process gg — > qq* which, at the tree level and in the massless 
limit, is described by the following matrix element squared: 

~ Nl-l' ■ ^^^^ 
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This does not have a unique rearrangement, not all such decompositions 
are positive, and many of them are not invariant under permutations of the 
external legs. 

As a further illustration of this second point, consider the process qq qq, 
which has two colour flow structures as shown in figure and whose matrix 
element squared is rearranged, according to as follows: 



Nc 



c 



+ 



+ 

2 



s^ + t' 

V? 

s 

' t 



c 



s 

ut 



(13) 
(14) 



This arrangement does not preserve the crossing symmetries (s ^ u) 
and {s ^ t) of the planar terms. It is in fact impossible to preserve the 
crossing symmetry while also preserving the pole structure, since the only 
dimensionless combination of s, t and u that is (s u) symmetric is su/t^ . 
Taking into account the pole structure, the nonplanar term /ut can not be 
expressed as a linear combination of su/t^, st/v^ and a constant. 

Although ( p!^ seems a natural rearrangement, it fails when Nc = 2 and 
the leg permutation (s ^ u) is made, corresponding to qq — > qq: 
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If we now set Nc ■ 
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2 in the second term above, we obtain 

u^ + t^ u _ {-t){u-t) 
s 

which is negative whenever {u — t) is negative. For Nc 
expression remains positive but this is accidental^. 
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tf + -t^ > 0. 



(15) 



(16) 

3 the above 
(17) 
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One method for obtaining a positive definite decomposition is to split the 
nonplanar term as follows: 



It 



U 



s-t 



1 1 
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Since both contributions are positive, the full terms are also positive. How- 
ever, the fact remains that the decomposition is not universal under permu- 
tations of the external legs. 

■^In general, when a full term is positive for a certain Nc, one can prove trivially that 
it is also positive for all higher Nc since the planar term must always be positive. 
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3 Proposed new procedure 



To solve the above problems, we propose the following procedure. Let the 
overall matrix element squared be given by: 

Wtot = I ) + (^-P-) = ^planar + (^-P-) • (19) 



where |A^Pj is the planar term for the i-th colour flow, and (n.p.) represents 
the nonplanar part. Then the i-th full term is given by: 



l-^Pfull. = l-^Ptot (20) 

planar 



\M\\ + =^=^{n.p.). (21) 

planar 



\M\^ 



In other words, we split the nonplanar part by the ratios of the planar terms. 

This is positive definite, as is obvious in equation (|20|) . This method 
resolves all three drawbacks of the MW method listed above, and carries 
some additional advantages: 

• the correct pole structure in each full term is automatically ensured — 
the larger the planar term the larger the full term. No extra pole is 
introduced since the sum of the planar terms is positive definite (even 
in unphysical regions) and furthermore can not approach zero faster 
than the nonplanar part; 

• related processes always have the same decomposition; 

• we do not even need a compact analytical expression for the matrix 
elements squared to carry out the colour decompositions. 

This last feature becomes useful for more complex n-body processes^], where 
the most efficient method for calculating matrix elements may be to utilise 



the helicity amplitude formalisms [13 



One drawback of this method is that the formulae may not always be 
aesthetically preferable to those obtained with the method of Marchesini 
and Webber, for example in the case ofQ qq — qq. However, this is not a 
problem in computer simulations, which are the only circumstances where 
the methods described in this paper are put to practice. 

^In a variant of this procedure is adopted as the default in the process e+e~ qqgg- 
^In fact, for aU other massless QCD processes the procedure described here yields results 
identical to those in [Q. 
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4 The 2 parton 2 sparton processes 

The gauge invariance and the simphcity of these processes greatly compactify 
the expressions for matrix elements. In particular, as in the QCD case, all 
supersymmetric processes of the form qq — ^ gg (where q and g refer here to 
quarks or squarks and gluons or gluinos, respectively) can be expressed as: 



(colour factor) x 



c;n.p. 



X |7W|2 



QED 



[22) 



provided that the sparticles exchanged in the t and the u 
masses m3(4) and m4(3) respectively for the process 1 + 2 - 



channels have 
^3 + 4. Here 



M4 



U 



-2pi ■ Pi and h = t 



rrin 



—2pi ■ p3. When rrit ^ ^3(4) or 
vtiu 7^ ?7i4(3), as is the case generally in qq ^ gg, there are correction terms 
proportional to the differences in the squared masses. 

The colour flows are identical to the ordinary QcD case shown in figures 



We present the spin and colour averaged squared matrix elements, gen- 
eral for any Nc and nondegenerate squark masses. These are divided by a 
statistical factor of two when the final state spartons are identical. 

Apart from the colour structures, the formulae match those of [|15[ when 
Nc = 3, and those of tree-level expressions in when left and right squark 
masses are taken to be degenerate. Stop and sbottom mixings are not con- 
sidered here since we are dealing with chirality independent interactions. 

In the formulae that follow, gs is the strong coupling, evaluated at some 
scale which is not determined at tree level (in the HERWIG Monte Carlo it is 
taken to be i/s). s is the effective centre-of-mass energy s^e = 'Stota;ia;2, t 
and u similarly. Since we are taking the initial state partons to be massless it 
follows that s+t+u = m\+m\ and s+t^+u^ = 0. Also ut—m\m\ = sp^ > 0, 
where Pt is the outgoing transverse momentum, is the gluino mass, 



are the i-ih. generation squark masses. Charge conjugate final states must 
also be included in simulations. 

gtc, 



\M\^{qiqi -^gg) 
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EC-. 

L,R 



L,R 



(23) 



with 



L,R 



2sp'^ 



(^1 - A^), + (tj - A\ - (.VAr^)^,^; 
(m4- A)(t3- A) 



ih - A)2 



(m4-A)2;, Nl.\ts-A Ui-Aj^ 



(24) 
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where A 



mi 



r denotes a summation over the left and right 



squarks. Terms marked with subscript t correspond to the colour flow (1 
3,3 4,4 2) and those with subscript u correspond to the colour flow 
(1 ^ 4, 4 — >• 3, 3 ^ 2), as in flgure|^. Although it is possible to distribute the 
nonplanar terms between the two colour flows using the 'intuitive' method, 
for the sake of consistency we advocate the use of equation (pl|). 



X 



,2 , ^2 , 4m?s2p2, 



ul + ti + 



sHl 



St 



1 



"4^3 



ut- 



(25) 



The colour flows are (1 3, 3 
3,3 2,2 1) for l/s'^ul, and 



4,4 



4,4 ^ 2,2 
[1 ^ 4,4 - 

(flgure §). There is no nonplanar term, as any cross term between two colour 
flows is equivalent to the square of the other. 



-> 1) for l/sHl (1 
2, 2 ^ 3, 3 ^ 1) for l/tlu'^ 



2 2 

-M4 - 2{ml - m^) ( 1 + ^ + ^ 



^, {uDs + {S')u - (ti/iVg)n.p. ^26) 
StsU4 

The colour flows are (2 — 1, 1 — 3, 3 — 4) for the s-channel term and 



(2— >-3,3^1,1^4) for the -u-channel term (flgure |^). Equation (21) leads 
to the decomposition of the nonplanar term t| = [{ultl)s + {sHl)u]/ {u1 + s^). 



2Nc ' 1 + 



The colour flows are (1 — > 4, 2 — * 3) for the /-channel term and (1 ^ 3, 2 — >■ 4) 
for the M-channel term (flgure 0). For the case i j, the colour flow is 
uniquely (1 — 4, 2 ^ 3), as in flgure |l|. 



\M\'^{q^qj q^^^q. 



3R,L/ 
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The colour flows are again (1 ^ 4, 2 ^ 3) for the t-channel term and (1 — > 
3, 2 — > 4) for the w-channel term. There is no nonplanar term. For the case 
z 7^ j, the colour flow is again uniquely (1 — > 4, 2 — > 3). 



2N, 



c 



{t - mlf 



2/N, 



c 



s{t — wA 



n.p. 



(29) 



The colour flows are (1 — > 2, 4 — 3) for the t-channel term and (1^3,4^2) 
for the s-channel term. For the case i ^ j, the colour flow is uniquely 
(1^2,4^3). 



\Mnm, - g.,,,4,J = 1^ ■ ^^-^ (30) 
The colour flow is uniquely (1 — > 2, 4 — ^ 3), regardless of i and j. 



|A^P(M.-W4,J = ^-^ (31) 



Here i ^ j. The colour flow is uniquely (1 3, 4 ^ 2). 



(32) 



{ul)t + {tl)u - {sVNlU. 



sHlul 



The colour flows are (4 ^ 2, 2 — > 1, 1 ^ 3) for the t-channel term and 
(4 ^1,1^2, 2^ 3) for the u-channel term (figure §). 



5 Conclusions 

We have discussed the techniques for calculating the 'full terms' in QcD and 
supersymmetric QcD processes which involve multiple colour flows. The ex- 
act distribution of nonplanar terms between the full terms is expected to 
make little significant difference to the prediction of soft radiation patterns, 
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but we argued that the conventional method used to achieve this is inade- 
quate. 

We have presented the formulae for 2 parton 2 sparton processes, to- 
gether with the colour flows associated with them. These arc incorporated in 
the Monte Carlo event generator HERWIG 6.1, to be released shortly. However, 
the study of the showering and hadronisation of supersymmetric particles has 
not been carried out in sufficient detail, and it is clear that further investi- 
gations are necessary. 
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Figure captions 

[1] Colour flow for processes of the type qiq'2 —>■ qsQ^- 

[2] Colour flows for processes of the type qiq2 — > g3?4. 

[3] Colour flows for processes of the type gi^2 9^9 a- 

[4] Colour flows for processes of the type (71(72 9z9i- The direction of 
colour flows (indicateci by the arrows) is arbitrary. 
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